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ABSTRACT 

Cold Dark Matter theory predicts that the Local Group should contain many more dwarf-sized objects 
than the observed number of dwarf galaxies — the so-called sub-structure problem. We investigate 
whether the suppression of star formation in these small objects due to the presence of a photoionizing 
background can resolve the problem. We make use of results from recent hydrodynamic simulations 
to build a recipe for the suppression of gas infall into semi-analytic galaxy formation models, and use 
these to predict the luminosity function of dwarf galaxies in the Local Group. In the models without 
photoionization "squelching" , we predict a large excess of faint dwarf galaxies compared with the observed 
number in the Local Group — thus, the usual recipe for supernovae feedback used in semi-analytic models 
does not solve the sub-structure problem on its own. When we include photoionization squelching, we 
find good agreement with the observations. We have neglected tidal destruction, which probably further 
reduces the number of dwarf galaxies. We conclude that photoionizing squelching easily solves the sub- 
structure problem. In fact, it is likely that once this effect is taken into account, models with reduced 
small-scale power (e.g. Warm Dark Matter) would underproduce dwarf galaxies. 

Subject headings: cosmology: dark matter — galaxies: formation — galaxies: dwarf 
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CZ2 Cold Dark Matter (CDM) models with scale-invariant 
Sitial fluctuation spectra seem to provide a remarkably 
"successful framework for explaining a broad range of ob- 
servations pertaining to structure on scales larger than ~ 1 
^^pc. Only recently have N-body techniques achieved the 
JrKmamic range necessary to resolve sub-structures within 
other bound, virialized structures. This new generation 
of very high resolution simulations has demonstrated that 
the presence of numerous sub-structures is an unavoidable 
consequence within the cosmological constant (A) domi- 
nated "ACDM" models that are now very widely used. 
Such simulations predict that dark matter halos with the 
mass of the halos of the Milky Way or M31 contain within 
their virial radius several hundred objects with internal 
velocities similar to those of observed Local Group dwarf 
galaxies (Klypin et al. 1999 (KKVP); Moore et al. 1999 
(M99)), whereas only about 40 such galaxies are known in 
the Local Group. This result has come to be known as 
the "sub-structure problem" , and is a direct consequence 
of the relatively large amount of power on small scales in 
"standard" ACDM models. 

A variety of solutions have been proposed which solve 
the problem by modifying some basic tenant of the usual 
CDM canon. These include reducing the small-scale power 
by appealing to inflationary models with broken scale in- 
variance (Kamionkowski & Liddle 2000), or changing the 
nature of the dark matter, e.g. Warm Dark Matter (Hogan 
& Dalcanton 2000; Bode et al. 2001), Self-interacting Dark 
Matter (Spergel & Steinhardt 2000), or Strongly Annihi- 



lating Dark Matter (Kaplinghat, Knox & Turner 2000). It 
appears, however, that each of these solutions may vio- 
late other constraints (Barkana, Haiman, & Ostriker 2001; 
Gnedin & Ostriker 2001; Hui 2001; Miralda-Escude 2000). 
Therefore, in this paper, we investigate a less exotic solu- 
tion. 

It has long been appreciated that the accretion and cool- 
ing of gas by low mass halos will be "squelched" in the 
presence of a strong photoionizing background. This pro- 
cess has been studied by many authors using analytic ap- 
proaches (Ikeuchi 1986; Rees 1986; Babul & Rees 1992; 
Efstathiou 1992; Shapiro, Giroux, & Babul 1994), 1-D hy- 
drodynamic simulations (Haiman, Thoul, & Loeb 1996; 
Thoul & Weinberg 1996) and full 3-D hydro simulations 
(Quinn, Katz, & Efstathiou 1996; Weinberg, Hernquist, 
& Katz 1997; Navarro & Steinmetz 1997; Gnedin 2000). 
These studies concur in the conclusion that after reioniza- 
tion, gas accretion is suppressed in halos with virial tem- 
perature less than some characteristic value. The Jeans 
Mass is an obvious candidate for this critical mass scale 
(e.g. Rees 1986). However, because dark matter halos 
in CDM models have complex and varied formation histo- 
ries, it is clear that some sort of time-averaged Jeans Mass 
is needed to determine the gas content of any given halo 
at any given time. The recent analysis of Gnedin (2000; 
GOO), based on hydrodynamic simulations, suggests that 
the 'filtering mass', which corresponds to the length scale 
over which baryonic perturbations are smoothed in linear 
theory, provides a good description of this time-averaged, 
effective Jeans Mass. The net result is that halos with virial 
velocities ,$, 30-50 km/s are able to collect a substantial 
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amount of gas only if they collapse before the universe is 
reionized. 

Bullock, Kravtsov, & Weinberg (2000; BKW00) inves- 
tigated whether this effect could resolve the sub-structure 
problem, using a model based on estimates of halo collapse 
redshifts obtained from an extension of the Press-Schechter 
(1974) approximation, combined with constraints from N- 
body results. They concluded that a combination of tidal 
destruction and photoionization squelching could very nat- 
urally solve the sub-structure problem in the local group. 
However, their model did not treat gas physics or star for- 
mation explicitly, instead relying on a few phenomenolog- 
ical parameters to determine which halos would host ob- 
servable galaxies. Chiu, Gnedin & Ostriker (2001) also 
concluded that photoionization squelching would reduce 
the slope of the galactic mass-function on the scale of 
dwarf galaxy halos, using scaling arguments based on the 
Press-Schechter approximation and results from the simu- 
lations of GOO. It remains worthwhile, however, to inves- 
tigate this process within the more detailed framework of 
semi-analytic galaxy formation modeling (e.g. Kauffmann, 
White, & Guiderdoni 1993; Somerville & Primack 1999; 
Cole et al. 2000). 

In this paper, we use the recipe provided by Gnedin 
(2000), based on hydrodynamic simulations, to incorpo- 
rate the effects of photoionization squelching into a semi- 
analytic model of galaxy formation. We examine the ef- 
fect on the luminosity function of dwarf galaxies in Milky 
Way-sized halos and re-assess whether this process can 
solve the sub-structure problem. Our models are based 
on Monte Carlo merger trees and include treatments of 
sub-halo merging, gas cooling, star formation, supernovae 
feedback, chemical enrichment, and stellar population syn- 
thesis. In all the models presented here, we adopt a 
"standard" ACDM cosmology with il m = 0.3, n A = 0.7, 
Hq = 70 km/s/Mpc, and og = 1. 

2. SEMI-ANALYTIC MODELS 

2.1. Merger Trees and Sub-structure 

Our models are based on Monte Carlo realizations 
of merger trees constructed using the extended Press- 
Schechter formalism and the method of Somerville & Ko- 
latt (1999). We track the decay of the orbits of sub-halos 
due to dynamical friction after being subsumed within a 
larger halo, and when a sub-halo reaches the center of the 
host halo, it is assumed to merge with the central object 
(for details see Somerville & Primack 1999; hereafter SP). 
We neglect tidal destruction of sub-halos, and also satellite- 
satellite merging. We would first like to check that this 
approximate treatment of the survival of substructure is in 
agreement with N-body simulations. Fig. fu shows the cu- 
mulative velocity function predicted by the semi-analytic 
models for sub-halos that reside within the virial radius 
of a host halo with a mass of 1O 12 M . This is com- 
pared with the velocity function of sub-halos within 200 
h _1 kpc, for halos of similar mass from the N-body simula- 
tions of M99 and from recent simulations produced using a 
new version of the Adaptive Refinement Tree (ART) code 
(Kravtsov, Klypin, & Khoklov 1997) with variable mass 



as 'the ART simulations', are similar to those presented in 
KKVP99, but have higher mass resolution. The results are 
binned in terms of Vm ax , the circular velocity at the peak 
of the rotation curve, which is the quantity that is most 
robust for sub-halos identified in the simulations. We use 
the results of Bullock et al. (2001) to calculate Knax for 
the halos in the semi-analytic merger trees as a function of 
their mass and collapse redshift. The semi-analytic trees 
have been run with the same mass resolution as the N-body 
simulations (the smallest halos have a mass of 10 /j _1 M q ). 




(km/s) 



Fig. 1. — The volume density of sub-halos, as a function of their 
maximum circular velocity Vmax. The thick solid line shows the re- 
sults from many realizations of a semi-analytic merger tree (normal- 
ized to the volume within the virial radius). Thin solid lines with 
squares show the results for two halos from the simulations of M99. 
Dashed lines with diamonds show the results from the ART simu- 
lations (see text), and the straight dashed line shows the power-law 
fit to the simulations presented in KKVP. The N-body results are 
normalized to the volume within 200 h _1 kpc. 



We see that the semi-analytic merger trees produce more 
of the small sub-halos (K nax ;$ 20 km/s) than the ART 
simulations, as found by BKW00. The simulations of M99 
have a factor of ~2-3 more of the smallest halos than the 
ART simulations, but this is likely to be mainly due to 
the tt m = 1 SCDM cosmology used for the M99 simula- 
tions. The recent ACDM simulations presented by Font & 
Navarro (2001) and Governato et al. (2001) both appear 
to have about 40 % less substructure than the SCDM sim- 
ulations of M99, in good agreement with the ART results. 
Some of the discrepancy could also be due to differences in 
the algorithm used to identify sub-halos, or to the slightly 
different masses of the host halos. 

BKW00 adopted a model for tidal destruction of sub- 
halos in order to bring the number of surviving sub-halos 
in the semi-analytic merger trees into agreement with that 
found in the ART simulations. While it is clear that tidal 
destruction of sub-halos does occur in the simulations (A. 
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dwarf galaxies does occur in the Local Group (we are wit- 
nessing it in progress in the case of the Sagittarius dwarf) 
it is still possible that dissipationless simulations overes- 
timate the effects of tidal destruction relative to the real 
Universe, where the presence of baryons might make galax- 
ies more robust to disruption. We opt to neglect the effects 
of tidal disruption in this work, for this reason and also be- 
cause we wish to find out if the effects of photoionization 
squelching alone can solve the sub-structure problem. 
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2.2. Photoionization Squelching 

The basic ingredients of the models used here, which 
include gas cooling, star formation, supernovae feedback, 
chemical evolution and stellar population synthesis, are de- 
scribed in detail in SP and Somerville, Primack & Faber 
(2001). We introduce one new ingredient, the squelching 
of gas infall by a photoionizing background. We assume 
that the Universe becomes reionized instantaneously at 
a redshift z- lon , which we treat as a free parameter. If a 
halo collapses before z- lon , it is allowed to capture a mass 
M g = /{,M v i r of gas which subsequently becomes available 
for cooling and star formation, where /t> is the universal 
baryon fraction and M v ; r is the virial mass of the halo. Af- 
ter reionization, the gas captured by a halo is given by the 
fitting function provided by Eqn. 7 of GOO: 



M g = ¥ 



f b M v 



0.26M C /M V 



(1) 



where Mc is a characteristic mass, which represents the 
mass of objects that on average retain 50% of their gas. 
We approximate Mc as the mass corresponding to a halo 
with virial velocity Vc — 50 km/s collapsing at redshift z 
(see the appendix of SP). This provides a good description 
of the redshift dependence of Mc found in the simulations 
of GOO: Mc increases from about 10 8 M Q at z — 8 to 
several 10 10 M© at z = 0. Then, in agreement with the re- 
sults shown in Fig. 2 of GOO and with previous results (e.g. 
Thoul & Weinberg 1996; Quinn et al. 1996), we find that 
halos with virial velocity V c ~ 50 km/s are able to cap- 
ture about half of the gas that they would accrete in the 
absence of a background, while halos with V c < 30 km/s 
are able to accrete very little gas, and halos with V c > 75 
km/s are nearly unaffected. 
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Fig. 2. — The luminosity function of the Local Group. Symbols 
with error bars show the observed local group luminosity function, 
with Poisson errors. The dotted line shows the results of the models 
with no feedback of any kind, and the solid line shows the results of 
the models with supernovae feedback but no photoionization squelch- 
ing. Short dashed, long dashed, and dot-dashed lines show the model 
results with squelching, for three different assumed epochs of reion- 
ization as indicated on the figure panel. The magnitudes of several 
familiar Local Group members are indicated. 

We normalize our models to reproduce the V-band lu- 
minosity of the central Milky Way or M31-like galaxies in 
our "Local Group" , which we represent as two equal mass 
halos each with M v - lr = 10 12 A/ Q . The luminosity function 
of dwarf satellite galaxies for an ensemble of such halos is 
shown in Fig. |[ The data points show the V-band lumi- 
nosity function of dwarf galaxies in the Local Group from 
the compilation of Irwin 2 (see also Mateo 1998). We show 
the model results without supernovae (SN) feedback or 
squelching, and for the models including supernovae feed- 
back according to the usual treatment (e.g. Kauffmann, 
White & Guiderdoni 1993) but without photoionization 
squelching. For the models without SN feedback, the fall- 
off at My ;$ —13 is due to the fact that gas within these ha- 
los cannot cool efficiently via atomic processes, and we have 
not included cooling by molecules, which we expect to be 
inefficient (Machacek, Bryan & Abel 2001; Haiman, Abel 
& Rees 2000). Including SN feedback pushes this "bend" 
to fainter magnitudes. The models including supernovae 
feedback only, at a level similar to previously published 
semi-analytic models, fit the luminosity function of satel- 
lite galaxies down to about My ~ —15 (a bit fainter than 
the SMC), but overproduce fainter objects by a consider- 
able factor — more than one order of magnitude for the 
faintest objects (such as Draco and Carina). 

We now consider the effect of photoionization squelching 
as discussed above. Results are shown for three different 
values of the reionization epoch z- mn , and are seen to be 
relatively insensitive to this parameter within a reasonable 
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range of values (note that observations constrain the epoch 
of reionization to 6 <; z ro ion <S 40, and most CDM-based 
modelling suggests that it lies in the range 6 <J z ro ion ^ 20; 
see the recent reviews by Shapiro (2001) and Barkana & 
Loeb (2001) and references therein). The large galaxies 
(the "Milky Way" and "M31") are virtually unaffected, 
the LMC-sized galaxies are slightly affected, and galaxies 
smaller than the SMC are significantly affected by the pres- 
ence of the background. The models now produce a rea- 
sonably good fit to the observed dwarf satellite luminosity 
function, especially given that tidal destruction may fur- 
ther reduce the number of satellites and that the observed 
LF may be incomplete for the faintest objects. 

3. DISCUSSION AND CONCLUSIONS 

We have compared the multiplicity function of sub- 
halos within a "Milky Way" sized halo predicted by semi- 
analytic Monte Carlo merger trees with that found in 
high-resolution N-body simulations. In agreement with 
BKW00, we found that semi-analytic merger trees predict 
more substructure than is found in the ART simulations. 
This is probably due to the tidal destruction of sub-halos 
in the N-body simulations, which we have chosen to ne- 
glect in the semi-analytic models. In continuing to assess 
the sub-structure issue, it will be important to determine 
the importance of tidal destruction of galaxies when the 
presence of baryons and the effects of baryonic contrac- 
tion on the halo are included. The simulations of M99 
have a factor of 2-3 more sub-halos than the ART sim- 
ulations, which is probably due to the SCDM cosmology 
used by M99. The physical origin of this cosmology depen- 
dence (now confirmed by several groups) is not yet clear 
and should be investigated. 

All previously published studies of which we are aware 
(e.g. M99; KKVP99; BKW00) have illustrated the sub- 
structure problem in terms of the circular velocity function 
of sub-halos and satellite galaxies. However, in order to do 
this the observed velocity dispersions of the dwarf satel- 
lite galaxies (which are predominantly spheroidal) must be 
converted to rotation velocities by making some assump- 
tions about the geometry and the potential. It is generally 
assumed that the orbits are isotropic and the potentials 
isothermal (M99); however, this assumption is uncertain 
and probably inaccurate at the 50% level, as the veloc- 
ity dispersions of many spheroids drop in the inner region 
where the dwarf velocity dispersions are measured (White 
2000) . In addition, the presence of baryons modifies the ro- 
tation curve, so that the peak velocity of the dark matter 
halo does not necessarily correspond with the peak velocity 
of the galaxy that forms within it. Moore (2001) and Font 
& Navarro (2001) attempt to account for these effects, and 
it appears unlikely that these corrections are large enough 
to make the sub-structure problem go away. Still, it is use- 
ful to present the comparison in a different form, in terms 
of the predicted and observed luminosity function of dwarf 
satellites. While many uncertain ingredients go into mod- 
cling galaxy luminosities, at least this quantity is directly 
observed and the selection effects are relatively well un- 
derstood. We find, in agreement with the previous studies 
based on velocity functions, that models without any kind 



satellite galaxies in the Local Group. 

We find that the standard implementation of supernovae 
feedback generally used in semi-analytic galaxy formation 
models cannot, by itself, solve the substructure problem. 
If we increase the efficiency of SN feedback to the point 
where enough dwarf galaxies are suppressed to solve the 
sub-structure problem, the feedback produces too great 
an effect on larger galaxies, introducing a curvature in 
the Tully-Fisher relation that is incompatible with obser- 
vations (see SP; also Moore 2001). Supernovae feedback 
might be made to work by introducing a more complicated 
mass dependence (rather than a power-law with a single 
slope), but this would have to be formulated in a purely 
ad-hoc way, and still might have difficulties with other ob- 
servational constraints. 

When we include the suppression of gas accretion by 
small halos after reionization in our semi-analytic model of 
galaxy formation, using a recipe based on state-of-the-art 
hydrodynamic simulations, we find that with no additional 
tuning, the model then produces very good agreement with 
the observed luminosity function of satellite galaxies in the 
Local Group. The properties of luminous galaxies in the 
models (which have been shown previously to agree well 
with observations) are nearly unaffected. The results are 
not very sensitive to the redshift of reionization, which we 
treat as a free parameter, although squelching is slightly 
less effective (leading to survival of more dwarf galaxies) 
when reionization is assumed to be later, because more 
galaxies are able to form before the squelching is turned 
on. 

It should also be noted that here we have considered only 
one possible form of photoionization feedback, the suppres- 
sion of gas collapse due to a uniform background produced 
by distant stars or quasars. We have neglected several fur- 
ther effects. Even gas that has managed to collapse might 
be "boiled" out of very small halos (Barkana & Loeb 1999; 
Haiman, Abel & Madau 2001). This effect is expected to 
be important mainly in the smallest of the halos harbor- 
ing observable dwarf galaxies (V c £ 20 km/s) . In addition, 
the presence of the background would reduce the efficiency 
of atomic cooling. Furthermore, the effect of turning on a 
quasar or starburst in close proximity to the dwarf halos (as 
for example, if our Galaxy or M31 underwent such episodes, 
which presumably they did at some point in the past, since 
both have bulges and probably black holes) might lead to a 
more dramatic form of photoionization feedback than the 
one we have considered here. Some combination of these 
effects might help to reduce the remaining excess of very 
faint (My ~ — 8) satellites in our models. 

For My ,$, — 10, the models with squelching come very 
close to producing the correct number of dwarf galaxies 
even neglecting the above processes and with no tidal de- 
struction. If our treatment of the efficiency of squelch- 
ing, supernovae feedback, and star formation is reasonably 
accurate, this suggests that tidal destruction may not be 
as efficient as dissipationless simulations (e.g. BKW00) 
would suggest — otherwise, our models would not pro- 
duce enough dwarf galaxies! Similarly, models with re- 
duced power on small scales (such as Warm Dark Matter) 
may actually suffer from a reverse sub-structure problem, 
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and fail to produce a sufficient number of dwarf satellite 
galaxies. 

The idea that we have explored here is very similar to 
that investigated by BKWOO. However, BKWOO did not 
explicitly model any gas processes (cooling, star forma- 
tion, SN feedback, etc) and so were forced to introduce 
a more simplified criteria to determine which sub-halos 
would harbor observable galaxies. While our models use 
simple recipes to represent the relevant physics, they do 
contain a self-consistent treatment of the many interwo- 
ven processes that determine the present day luminosity of 
a dwarf galaxy. We therefore consider it encouraging that 
our investigation has produced similar conclusions to those 
of BKWOO. 

While there have been suggestions that the photoion- 
ization phenomenon could be responsible for producing 
Lyman-limit QSO absorption systems at z > 3 (Abel & 
Mo 1998) and bursting dwarf galaxies at z ~ 1 (Babul & 
Rees 1992), and the effects of photoionization feedback on 
the epoch of reionization itself have been explored (Haiman 
& Loeb 1997; Ciardi et al. 2001), there have been rela- 
tively few investigations connecting this physics with the 
properties of nearby dwarf galaxies, perhaps because most 
semi-analytic galaxy formation models have not included 



this effect to date. With this mechanism in place, we are 
now in a position to investigate many other properties of 
dwarf galaxies within this framework. 

The squelching scenario has a number of interesting im- 
plications. It suggests that there are many barren halos 
floating around in the Local Group, as well as many very 
faint galaxies that lie below current detection limits. There 
has been some concern that these halos might heat the disk 
of our Galaxy beyond allowed observational limits, but re- 
cent studies (Velazquez & White 1999; Font & Navarro 
2001) suggest that the halos in question are too small to do 
much damage. However, the prospects of detecting these 
halos via gravitational lensing or by direct detection if the 
Earth should happen to pass through one of them (Moore 
et al. 2001) remain intriguing possibilities. 

ACKNOWLEDGMENTS 

We thank B. Moore, A. Klypin, and A. Kravtsov for mak- 
ing the results of their simulations available to us electronically, 
and for useful discussions. We also thank R. Barkana for useful 
discussions, and T. Abel, J. Bullock, and A. Kravtsov for help- 
ful comments on the manuscript. RSS is supported by a rolling 
grant from PPARC. 



REFERENCES 



Abel, T. & Mo, H.J. 1998, ApJ, 494L, 151 

Babul, A., & Rees, M.J. 1992, MNRAS, 255, 346 . 

Barkana, R., Haiman, Z. & Ostriker, J.P. 2001, kstro-ph/01023a4 
Barkana, R._& Loeb, A. 2001, Physics Reports, in press, isto 

ph/00104«S 

Bode, P., Ostriker, J.P., & Turok, N. 2001, ApJ, 556, 93 
Bullock, J.S., Kolatt, T.S., Sigad, Y., Somerville, R.S., Kravtsov, 

A.V., Klypin, A., Primack, J.R., & Dekel, A. 2001, MNRAS, 321, 

559 
Bullock, J.S., Kravtsov, A.V, Weinberg, D.H., 2000, ApJ, 548, 33 

(BKWOO) 
Cole, S., Lacey, C.G., Baugh, CM., & Frenk, OS. 2000, MNRAS. 

319, 168 
Chiu, W., Gnedin, N.Y., & Ostriker, J.P. 2001 
Ciardi, B. Ferrara, A., Governato, F., Jenkins 

2000, 314, 11 



istro-ph/0103359| 
A. 2000, MNRA S, 



Efstathiou, G. 1992, MNRAS, 256, 



-4_L 



astro-ph/0106268 



Font, A.S., & Navarro, J.F. 2001, 

Gnedin, N.Y. 2000, ApJ, 542, 535\v»w^ . 

Gnedin, O.Y. & Ostriker, J.P. 2001, astro-ph/0010436| 

Haiman, Z., Abel, T., & Rees, M.J. 2000, 534, 11 

Abel, T., & Madau, P. 2001, ApJ, 551, 599 
& Loeb, A. 1997, ApJ, 483, 21 
Thoul, A.A., & Loeb, A. 1996, ApJ, 464, 523 
Dalcanton, J.J. 2000, Phys. Rev. D62, 063511 

Hui, L. 2001, Phys. Rev. Lett. 86, 3467 

Ikeuchi, S. 1986, AP&SS, 118, 509 



Haiman, Z.. 
Haiman, Z.. 
Haiman, Z.. 
Hogan, C.J. 



astro-ph/0005210 



Kaplinghat, M., Knox, L., & Turner, M.S. 2000, 
Kamionkowski, M. & Liddle, A.R. 2000, Phys. Rev. Lett. 84, 4525 
Kauffmann, G., White, S.D.M., & Guiderdoni, B., 1993, MNRAS, 
264, 201 



Khokhlov, A.M. 1997, ApJS, 111, 



Klypin, A.A., Kravtsov, A.V., Valenzuela, O., & Prada, F. 1999, 

ApJ, 522, 82 (KKVP) 
Kravtsov, A.V., Klypin, A. A 

73 
Machacek, M.E., Bryan, G.L., & Abel, T. 2001, 548, 509 

Mateo, M.L. 1998, ARA&.A , 36, i35 . 

Miralda-Escude, J., 2000, kstro-ph/0002050 

Moore, B., Calcajieo-Rolc han; C , SB3S1; ', Quian , T., Lake, G., 

Ghigna, S., Governato, F. 2001, 3,stro-ph/0106271 
Moore, B., 2001, Plenary talk, 20tn lexas symposium, rds J C. 

Wheeler & H. Martel, AIP Conference Series, kstro-ph/0103100 
Moore, B., Ghigna, F., Governato, F., Lake, G., Stadel, J., & lozzi, 

P. 1999, ApJ, 524, L19 (M99) 
Navarro, J., & Steinmetz, M. 1997, ApJ, 478, 13 
Press, W.H., & Schechter, P. 1974, ApJ, 187, 425 
Quinn, T., Katz, N, & Efstathiou, G. 1996, ApJ, 278, 49 
Rees, M.J., 1986, MNRAS, 218, 25 
Shapiro, P. 2001, Plenary Session, 20th Texas Sy mposium, eds . J ,C. 

Wheeler & H. Martel, AIP Conference Series, astro-ph/0104315 
Shapiro, P.R., Giroux, M.L., & Babul, A. 1994, ApJ, 427, 25 
Spergel, D.N., & Steinhardt, P.J. 2000, Phys.Rev.Lett., 84, 3760 
Somerville, R.S. & Kolatt, T.S. 1999, MNRAS, 305, 1 
Somerville, R.S. & Primack, J.R. 1999, MNRAS, 310, 1087 
Somerville, R.S., Primack, J.R., & Faber, S.M., 2001, MNRAS, 320, 

504 
Thoul, A.A. & Weinberg, D.H. 1996, ApJ, 465, 608 
Velazquez, H., & White, S.D.M. 1999, MNRAS, 304, 254 
Weinberg, D.H., Hernquist, L., & Katz, N. 1997, ApJ, 477, 8 
White, S.D.M. 2000, ITP conference, urkonline. itp.ucsb.edu/online/galaxy_c00 



